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Abstract 

Momenta and masses of heavy projectile fragments (Z > 8), produced in 
collisions of 197 Au with C, Al, Cu and Pb targets at E/A = 600 MeV, were 
determined with the ALADIN magnetic spectrometer at SIS. Using these in- 
formations, an analysis of kinematic correlations between the two and three 
heaviest projectile fragments in their rest frame was performed. The sensi- 
tivity of these correlations to the conditions at breakup was verified within 
the schematic SOS-model. For a quantitative investigation, the data were 
compared to calculations with statistical multifragmentation models and to 
classical three-body calculations. 

With classical trajectory calculations, where the charges and masses of 
the fragments are taken from a Monte Carlo sampling of the experimental 
events, the dynamical observables can be reproduced. The deduced breakup 
parameters, however, differ considerably from those assumed in the statistical 
multifragmentation models which describe the charge correlations. If, on the 
other hand, the analysis of kinematic and charge correlations is performed for 
events with two and three heavy fragments produced by statistical multifrag- 
mentation codes, a good agreement with the data is found with the exception 
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that the fluctuation widths of the intrinsic fragment energies are significantly 
underestimated. A new version of the multifragmentation code MCFRAG 
was therefore used to investigate the potential role of angular momentum at 
the breakup stage. If a mean angular momentum of 0.75/i/nucleon is added 
to the system, the energy fluctuations can be reproduced, but at the same 
time the charge partitions are modified and deviate from the data. 

PACS number(s): 25.70.Mn, 25.70.Pq, 25.75.Ld, 25.75.-q 
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I. INTRODUCTION 



In several experiments with the ALADIN spectrometer, the decay of excited projectile 
spectator matter at beam energies between 400 and 1000 MeV per nucleon was studied . 
In these collisions, energy depositions are reached which cover the range from particle evapo- 
ration to multi-fragment emission and further to the total disassembly of the nuclear matter, 
the so-called 'rise and fall of multifragment emission' Q. The most prominent feature of 
the multi-fragment decay is the universality that is obeyed by the fragment multiplicities 
and the fragment charge correlations. These observables are invariant with respect to the 
entrance channel - i.e. independent of the beam energy and the target - if plotted as a 
function of Zb oun d, where Zbound is the sum of the atomic numbers Zi of all projectile frag- 
ments with Zi > 2. For different projectiles, the dependence of the fragment multiplicity on 
Zbound follows a linear scaling law. These observations indicate that compressional effects 
are only of minor importance. In contrast to central collisions at lower energies, where large 
radial flow effects are observed, the quantitative interpretation of kinematic observables is 
therefore simplified. 

More important, these characteristics are an indication that chemical equilibrium is 
attained prior to the fragmentation stages of the reaction. In fact, statistical models were 
found to be quite successful in describing the experimental fragment yields and charge 
correlations, if the breakup of an expanded system was assumed |fH^ H^?l - In addition, 
the temperature of the excited matter, extracted from double ratios of isotope yields, is 
reproduced. On the other hand, the kinetic energy spectra of particles and fragments are 
not equally well described within the statistical picture. The energy spectra of light charged 
particles (A < 4) can be explained by a thermal emission of the fragments, but their slopes 
correspond to temperatures approximately three times larger than those extracted from 
isotope ratios While this may be an indication for pre-breakup emission if remains to 
be understood whether the kinetic energies of intermediate mass fragments (3 < Z < 30) 
are consistent with the statistical approach. 

The dynamics of the multifragmentation process has therefore to be studied. It is well 
known that kinematic correlations, which are governed by the long range Coulomb repulsion, 
are sensitive to the disintegration process. Previous studies concentrated mostly on the two- 
fragment velocity correlation functions ]T4|-f2"2]| . Only few attempts were made to analyze 
higher order correlations. However, these studies were done either for heavy fragments at 
much lower beam energies |23|-p7f or for light charged particles only ||28|| . In this paper, 
the results of a kinematic analysis of the fragmentation process of the projectile spectator 
are presented. Heavy projectile fragments produced in peripheral Au induced collisions at 
E/A = 600 MeV are studied without the influence of energy thresholds of the detectors. 
Moreover, the analysis is performed in the center of mass frame of the fragments, thus 
reducing the influence of directed collective motion of the emitting source. On the other 
hand, a limit of Z > 8 is imposed for kinematic observables by the lower detection threshold 
of the TP-MUSIC II tracking detector. The analysis is therefore restricted to the excitation 
energy range characterized by increasing fragment multiplicities. 



II. THE EXPERIMENT 
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A. Experimental setup 



The experiment was performed with the ALADIN forward spectrometer at the heavy-ion 
synchrotron SIS of the GSI Darmstadt, using a gold beam with an energy of 600 MeV per 
nucleon and a typical intensity of 2000 beam particles during a 500 ms spill. A schematic 
view of the experimental setup in the bending plane of the ALADIN magnet is shown in figure 
[l]. The incoming beam entered the apparatus from the left and first hit the beam counters, 
where for each beam particle the position in a plane perpendicular to the beam direction and 
the arrival time were measured with resolutions 5 X ~ 5 y ~ 0.5 mm FWHM and S t =100 ps 
FWHM, respectively. One meter downstream, the reaction target was positioned. Targets of 
C, Al, Cu and Pb with a thickness between 200 and 700 mg/cm 2 were used, corresponding 
to an interaction probability of up to 3%. Light charged particles from the mid-rapidity 
zone of the reaction were detected by a Si-CsI-array which was placed at angles between 
7° and 40° with a solid angle coverage of approximately 30% in this angular range ( 50% 
between 7° and 25°, 15% between 25° and 40° ). Fragments from the decay of the projectile 
spectator, emitted into a cone of approximately 5° around the direction of the incident 
beam, entered the magnetic field of the magnet. The magnet was operated at a bending 
power of 1.4 Tm which corresponded to a deflection of 7.2° for fragments with beam rigidity. 
The particles were detected in the TOF-wall, which was positioned 6 m behind the target. 
The time-of-flight of light particles with respect to the beam counter was measured with a 
resolution of 300 ps FWHM and with a resolution of 140 ps FWHM for particles with a charge 
of 15 and above. The TOF-wall provided the charge of all detected particles with single 
element resolution for charges up to eight. Charged particles with charges of eight and above 
were simultaneously identified and tracked by a time-projecting multiple-sampling ionization 
chamber TP-MUSIC II ( see section |II B| ), which was positioned outside the magnetic field 
between the magnet and the TOF-wall. To minimize the influence of scattering, of energy 
loss and of secondary nuclear reactions of the fragments after their production in the target, 
the spectrometer up to an entrance window in front of the ionization chamber was operated 
in vacuum. The components of the apparatus with the exception of the MUSIC detector 
have already been described in M. 



B. The MUSIC detector 

The TP-MUSIC detector is a time-projection multiple-sampling ionization chamber. If 
a charged particle passes through its active volume, an ionization track containing positive 
ions - which will drift to the cathode - and free electrons - which will move in the direction 
of the anodes - is produced. Due to the homogeneous electric field, the drift velocity of the 
electrons towards the anodes is independent of the position within the gas volume. Therefore, 
the distance of the primary particle track from the anode is proportional to the time the 



center of the electron cloud needs to reach the anode. The version TP-MUSIC II 291 which 
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was used in this experiment is shown in figure |2| []. It consists of three active volumes with 
the drift field in adjacent sections perpendicular to each other, two for the measurement 
of the horizontal and one for that of the vertical position and angle of the particle track. 
Each field cage has an active area of 100 cm (horizontal) times 60 cm (vertical) and a 
length of 50 cm. The horizontal field cages are both divided into two halves with a vertical 
cathode plane in the middle of the detector to reduce the maximum drift length and the 
high voltages necessary to provide the drift field. The chambers were operated at a high 
voltage of 150 V/cm, i.e. 7.5 kV for the horizontal and 9 kV for the vertical field cages, P10 
( 90% argon, 10% methane ) at a pressure of 800 mbar served as counting gas. To allow 
multiple sampling of the particle signals, each anode is subdivided into 16 stripes with a 
width of 3 cm each. 

The anode signals were recorded using flash ADCs with a sampling rate of 16 MHz. 
Together with a drift velocity of the electrons of approximately 5.3 cm//xs this corresponds 
to amplitude measurements at a step size of 3 mm in the direction of the drift. Since the 
drift time of the electron cloud is measured by each of the 16 anodes of a field cage i.e. at 
16 points along the beam direction (z-direction), the complete track information both in x- 
( field cage 1 and 3 ) and y-direction ( field cage 2 ) of the primary charged particle inside 
the MUSIC volume is available. The detector is operated outside the magnetic field volume 
of the ALADIN magnet, therefore the ionization track through the MUSIC gas is a straight 
line which is obtained by fitting the 16 track positions by three straight lines - one in each 
field cage. 

The position resolution has been estimated using the fact that the horizontal component 
of a track is determined with two separate field cages. The intersections of the measured track 
segments from the first and the third field cage with a virtual reference plane, positioned 
in the center of the vertical field cage and perpendicular to the z-direction, are calculated. 
The distance between these two points of intersection is a measure of the overall position 
and angle resolution of the detector. Its distribution is a gaussian with a width of 2.4 mm 
FWHM for particles with a charge of 20 and above which increases to approximately 12 mm 
at the detection threshold of Z = 8. These values are of the same order of magnitude as the 
effect of small angle scattering of the fragments in the counting gas of the MUSIC. 

The amplitude of the primary signal produced by a particular fragment is proportional 
to q 2 (3 2 , where (3 is the velocity of the particle and q is its charge state. Fragments from 
the decay of the projectile spectator are moving approximately with beam velocity. In this 
case, all fragments with nuclear charges up to 50 are fully stripped after passing through the 
target matter. They remain fully stripped in the detector gas, the primary signal is therefore 
proportional to the square of the nuclear charge of the particle. For particles with nuclear 
charges between 50 and 79, the mean charge exchange length in the MUSIC gas ( 8 cm and 
30 cm for Z = 50 and 79, respectively ) is small compared to the path length of the particle 
within the MUSIC detector. They reach their equilibrium charge state within the detector 
volume and the primary signal is proportional to the square of the effective charge. 

The amplitude of the primary signal decreases due to diffusion broadening ( proportional 



x The data published in were taken with the version III of the TP-MUSIC and a larger 

TOF-wall. 
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to the drift distance ) and due to impurities of the counting gas ( proportional to the square 
of the drift distance ). The amplitude measured at the anode is therefore dependent on the 
drift distance of the electron cloud. To determine the position correction the incident beam 
- i.e. particles with known Z and j3 - is swept across the field cages by varying the field of 
the magnet. In addition, the signals are corrected for the deviations from the beam velocity. 
This is essential for the charge resolution of binary fission fragments which have the widest 
distribution of laboratory velocities of all heavy fragments (Z > 8) from the decay of the 
projectile spectator. A charge resolution of 0.5 charge units FWHM is reached. This is 
demonstrated in figure |3] where a charge spectrum of the MUSIC-detector is shown. Since 
both the differences in pulse height for two neighboring charges and their fluctuations are 
proportional to Z, the charge resolution is independent of the charge of the fragment. The 
lower threshold for particle identification reached in this experiment is Z = 8. 



C. Momentum and mass reconstruction 

From the tracks of the charged particles measured behind the ALADIN magnet, the 
rigidity vector can be determined if the magnetic field is known. It was decided to fit 
the particle properties as a function of the measured track parameters rather than using 
a backtracing method because the latter is more time consuming at the analysis stage. 
For particles with a rigidity vector R, trajectories starting at the target position z ta rg and 
coordinates (x tar g, Vtarg) within the beam spot are calculated using the routines provided 
by the program package GEANT |j32| . The starting conditions are chosen from a five- 
dimensional grid with equidistant spacing for the variables x targ , Vtarg, R x /R and 
Ry/R. For a given magnetic field strength of the ALADIN magnet, the intersection (x mus i C , 
Vmusic) of each track with the reference plane of the MUSIC detector and its angle (m x , m y ) 
relative to this plane as well as the path length to this point are determined. The bending 
plane of the magnet is the horizontal x-z-plane, i.e. the main component of the magnetic 
field points to the direction of the y-axis, although the fringe fields can not be neglected, 
especially if the full geometric acceptance is used. 

Since a large range in iV/Z-ratios (0.7 - 1.5) and emission angles has to be covered, only 
40% of the grid points correspond to trajectories which reach the reference plane behind 
the magnet, all others end at the wall of the magnet chamber where they are lost. For the 
successful tracks, the three components of the rigidity vector together with the path length 
are fitted as the product of one-dimensional functions of five variables: the position (x mus i C , 
Vmusic), the angle m x and the target position (x targ , Vtarg)- The fit is done by means of an 
expansion in series of Chebychev polynomials for each variable. For a magnet which has 
virtually a dipole field, the most relevant terms are linear in x mus i C and m x for 1/R, R x /R 
and the path length, and linear in y mus i C for R y /R, but for a accuracy of the momentum 
reconstruction on the percent level, higher order terms can not be ignored. Under the 
assumption of an expansion up to third order, approximately 1000 individual contributions 
have to be calculated, which is not feasible. However, a particular term can be estimated 
by the size of the related expansion coefficient, since Chebychev polynomials are orthogonal 
within the interval from -1 to 1, and at the same time all their minima and maxima within 
this interval have the values -1 and 1, respectively. ( Strictly mathematically speaking, this 
is not correct. Among other conditions, the orthogonality relations can only be used if the 
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full parameter space is covered. This is not the case, since not all of the tracks reach the 
reference plane. ) In a second step, small terms are gradually suppressed until the x 2 of 
the fit has increased by 10%, thus reducing the total number from between 400 to 1000 in 
the first step ( depending on the highest order taken into account ) to 25 to 40 ( depending 
on the variable ). The fitting procedure is then repeated using only the remaining relevant 
terms which leads to slightly different expansion coefficients in the final result. 

Once this fitting procedure has been performed for each setting of the magnetic field used 
during the experiment, the reconstruction of the rigidity vector and of the path length is 
reduced to the evaluation of a set of polynomials. If a reasonable qualitiy of the reconstruc- 
tion can be archieved, this is a justification for the somewhat heuristic method to select the 
relevant contributions. The accuracy can easily be determined by calculating tracks with 
random start values - i.e. with starting coordinates at the target and for the rigidities not 
identical with the starting parameters used for the fitting procedure. The reconstruction 
is done for these tracks by evaluating the fit functions, and the input values are compared 
to the reconstructed ones. The mean deviations for rigidity and path length are a measure 
of the uncertainty caused by the reconstruction method itself. Clearly, the size of these 
deviations is on the one hand dependent on the mesh size of the grid of start values and 
on the other hand on the choice of the highest order taken into account for the expansion 
in Chebychev polynomials. Both quantities were optimized until the internal accuracy for 
all variables was better than 0.1% FWHM within the chosen range of rigidities between 1.2 
and 3.6 GeV/c. The final set of coefficients was obtained by fitting ^12000 tracks with a 
maximum order of 4 for each polynomial and a maximum of 6 for the sum of the orders 
within a term. 

A very similar procedure as described above can be used to estimate the expected errors 
due to the experimental resolution of the two position detectors in front and behind the 
magnet. A random offset of the order of the experimental uncertainties is added to the 
positions and slopes prior to the evaluation of the polynomials. Afterwards, the difference 
between the reconstructed values with and without random offsets is calculated. The mean 
value of these deviations is the resolution expected due to the experimental uncertainties. 
It was found that both an uncertainty of 0.7 mrad and of 3 mm produce an error in the 
rigidity of 1%. With the time and position resolutions given in the previous section, rigidity 
resolutions of approximately 1.2% and 3% can be expected for beam particles and medium 
heavy fragments with charge ^12, respectively. The quality of the rigidity reconstruction 
can be demonstrated by the rigidity distribution of beam particles passing through the 
apparatus without any nuclear interaction. Within all targets used in this experiment, gold 
projectiles reach their equilibrium charge state, providing particles with identical momenta 
and charge states 77 + , 78 + , 79 + , i.e. with rigidities which differ by 1.3% per charge state. 
With the carbon target, the influence of angular straggling within the target is small and 
negligible compared to the experimental errors due to the position resolution. In figure |], 
the rigidity distribution of beam particles after passing through the carbon target is plotted 
versus their x-position in the MUSIC reference plane. In this representation of the data, the 
three charge states ( equilibrium charge state distribution of 600 MeV/u gold in carbon: 
59% of the projectiles are fully stripped, 35% have a charge state of 78 + and 6% of 77 + ) 
are clearly visible, i.e. the rigidity resolution for heavy nuclei is approximately 1.3% FWHM 
which is in agreement with values expected from the resolutions of the individual detectors. 
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Using the reconstructed values for the rigidity and path length, the charge of the particle 
measured by the MUSIC detector and the time of flight given by the TOF-wall, the velocity 
and the momentum vector can be calculated for each charged particle detected both in the 
MUSIC and the TOF-wall, i.e. for particles with a charge of eight and above. The knowledge 
of velocity and momentum allows the calculation of the particle's mass. 

In figure ||, the mass spectra for the reactions Au+Al and Au+Cu are shown. Single mass 
resolution for charges up to 12 is obtained, corresponding to a mass resolution AA/A of ap- 
proximately 4.0% FWHM for light fragments. The dominant contribution to the uncertainty 
of the mass measurement 

AA j(AR\ 2 ( ' ATOF \ 2 

Till) + l 7 -^tof) ■ (1) 

is caused by the mass dependent error of the time measurement which is amplified by the 
factor 7 2 ( 7 2 =2.6 for 600 MeV/nucleon ). From this, a rigidity resolution of 2.4% FWHM 
can be deduced for light fragments. 



III. DATA 

The breakup dynamics of multifragmenting spectator matter will be reflected in the 
momenta of the fragments produced. Especially observables combining the kinematic infor- 
mation of two or more particles - e.g. relative velocities - are governed by the long-range 
Coulomb repulsion and are therefore sensitive to time scales of the decay and spatial prop- 
erties of the decaying source. Clearly, the breakup pattern will change with increasing 
excitation energy transferred to the spectator matter. From the analysis of reactions of gold 
projectiles with different targets and beam energies between 400 and 1000 MeV/nucleon it 
is well established |f[|,P|,|3"U] that the quantity Z bound - defined as the sum of the charges of all 



particles with charge two and above, which are emitted from the projectile spectator and 
detected in the TOF-wall - reflects directly the size of the spectator as well as the excita- 
tion energy transferred to the spectator nucleus. It was furthermore shown that the mean 
number of fragments produced in a reaction as well as other observables characterizing the 
populated partition space were independent of the target used, if they were investigated as 
a function of the quantity Zbound- Z boU nd is therefore used as a sorting parameter describing 
the violence of the reaction. 

As was discussed earlier, momenta and masses could only be reconstructed for particles 
with a charge Z of eight and above. It will be shown in the next section that events 
with two and more large fragments with Z > 8 cover the Z bound range from 30 to 70. The 
maximum mean number of intermediate mass fragments - defined as fragments with charges 
between 3 and 30 - is observed for a Z bound value of approximately 40. The dataset available 
covers therefore the range from peripheral collisions up to the region of maximum fragment 
production. 



A. Characterization of two- and three-particle events 

To show the characteristics of the event classes with two and three heavy particles with 
charge Z > 8, their reaction cross sections da / dZ bound are plotted in the upper panel of 



8 



figure Hfor the four different targets as a function of Zbound- In the following, events with two 
(three) fragments with charge Z > 8 are called binary (ternary). Binary events attributed 
to binary fission were excluded by the condition that either the lighter fragment is of charge 
below 20 or the sum of the two charges is smaller than 60. In the Zi^-plane, this region is 
well separated from the region of binary fission |34| . For comparison, the inclusive reaction 
cross sections - i.e. without conditions on fragment multiplicity and charge - are also shown: 
The binary and ternary events represent approximately 10% and 1% of the nuclear reaction 
cross section, respectively. In order to demonstrate that binary fission events as defined 
above populate an impact parameter region different from that of binary events without 
fission, the cross section for binary fission in the reaction Au+C is included. These events 
are obviously produced in very peripheral reactions. It had been shown earlier f|| that 
multifragment events evolve - with decreasing Zbound - from events with one heavy residue 
in the exit channel of the reaction and not from binary fission events. 

The sum of the charges of the two and three fragments is plotted in the lower panels. 
In ternary events, typically 80% of Z bound is contained in the charges of the three heavy 
fragments with average charges and masses of (Zi) = 22, 13, 10 and (A4) = 48, 29, 20, 
i=l,2,3. In binary events, the sum of the charges of the heavy fragments accounts on 
average for 75% of Zbound with a clear minimum at Zf, oim d=40, where the maximum mean 
number of intermediate mass fragments is observed. The average charges and masses for 
this event class are (Zj) = 26, 13 and (Ai) = 57, 27, i=l,2. 

It will now be demonstrated that these two event classes are representative subsets of 
the experimental data, i.e. that for a given Zbound value no evidence for a strong dependence 
on the number of heavy fragments is found. This means that other quantities defining an 
event do not show a close correlation between their mean values and the multiplicity of 
the heavy fragments if analyzed according to Z bound . Evidently, only observables can be 
used for this investigation which are not dominated by autocorrelations. The multiplicity 
of intermediate mass fragments for instance contains the number of all heavy fragments 
with a charge smaller than 30. The mean multiplicity of IMFs is therefore influenced by 
the selection criterion and will be significantly different for events with different numbers of 
heavy fragments in the exit channel. 

The mean number (Mi p ) of light particles from the mid-rapididy zone of the reaction 
which were detected in the hodoscope is a quantity which is certainly dependent on the 
violence of the reaction but independent of the specific decay channels of the excited pro- 
jectile spectator. In figure [7|, the inclusive distributions of (Mj p ) versus Z hound for the four 
targets are shown together with the distributions for events with two and three heavy par- 
ticles. In agreement with the participant-spectator-model, the size of the interaction zone 
- represented by the mean number of light particles - increases with decreasing size of the 
projectile spectator. The distributions are independent of the multiplicity of heavy pro- 
jectile fragments with the exception of the most peripheral reactions ( Z bound > 65 ). In 
this range of largest impact parameters the inclusive data are dominated by spallation and 
not by multifragment at ion events. There, the restriction to events with two or three heavy 
particles in the exit channel is synonymous with the selection of events with higher mean 
energy. 

The transversal deflection of the decaying projectile spectator is another quantity which 
is not influenced by autocorrelations with regard to the decay pattern. Since in events with 
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two or three heavy particles in the exit channel the heavy particles contain typically 75-80% 
of Zb oun d, the center of mass of these particles is in good approximation the center of mass 
of the decaying system. Thus, the transversal velocity 



of the center of mass of the two or three particles with respect to the beam frame was cal- 
culated. In figure ||, the mean values of this velocity as a function of Z bound are compared 
for events with two and three heavy fragments in the exit channel. In agreement with inclu- 
sive measurements at 400 MeV/nucleon P5] , jStrans increases monotonously with decreasing 
Zbound and establishes the transversal deflection of the projectile spectator and therefore the 
transversal momentum transfer (bounce) as a measure of the deposition of excitation energy 
into the spectator matter. Pure Coulomb interaction during a grazing collision would lead 
to very small values for the bounce between 5- 10~ 4 c and 4- 10 _3 c for C and Pb, respectively. 
But due to the trigger condition demanding at least one light particle detected in the ho- 
doscope and therefore a nuclear reaction, the bounce does not vanish for Zbound = 80. The 
increasing Coulomb repulsion with increasing charge of the target nucleus is nevertheless 
reflected in the small target dependence. Within the experimental errors, the transversal 
velocity at a given Z bound is independent of the two decay patterns studied. 

The two quantities (Mi p ) and Ptrans describe properties related to the initial reaction 
phase - the size of the fireball and the excitation energy transferred to the spectator matter. 
The fact that these quantities are independent of a specific choice of the multiplicity of heavy 
fragments demonstrates that a restriction to the subset of events, defined by the detection 
threshold of the MUSIC detector, does not select a non-typical sample of the produced 
projectile spectators. 



From the measured momenta of the heavy fragments the intrinsic momenta p C m{i) and 
velocities v cm (i) in the center of mass frame (CM- frame) of the binary or ternary heavy 
fragment system werde determined. This has primarily the advantage of eliminating the 
projectile velocity from the analysis. Furthermore, it reduces the influence of directed col- 
lective motion on the momenta of the particles. This is especially important if the data are 
to be compared to calculations with models which do not include linear collective motions. 
By construction, these momenta are collinear in the case of two and coplanar in the case 
of three particle events. For the further analysis, a new coordinate system has been chosen 
such that for each event the momentum vectors lie in the same plane - the xy-plane - and 
that the direction of the heaviest particle coincides with the x-axis. This eliminates the three 
Euler angles which describe the spacial orientation of the momenta relative to the beam axis. 
The kinematics of the two and three heavy fragments is thus reduced to one ( p x (l) ) and 
three ( Px(l), Py(2), Px(2) — p x (3) ) parameters, respectively. The relative kinematics of the 
fragments can thus uniquely be expressed in terms of one and three independent quantities 
which, for the analysis presented in this paper, are chosen as follows: (i) the total kinetic 
energy E 3 of the fragments in the CM-frame, (ii) the reduced relative velocity tVed(2, 3), and 
(iii) a quantity which describes the event shape in velocity space. In the case of only 
two heavy particles, the kinetic energy E 2 alone is sufficient to describe the decay dynamics. 




(2) 



B. Two- and three-particle observables 
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The sum E3 of the kinetic energies of the three particles is calculated in their CM-frame 



E, 



3 

3 = E 



Pc 



\ 2 • m ■ Ai 



(3) 



where m =931.5 MeV/c 2 is the atomic mass unit and A4 the mass number of the fragment 
i. The kinetic energy of the particles is dominated by the Coulomb interaction which itself 
is strongly dependent on the charges involved. The mean value (E 3 ) is therefore studied 
together with the standard deviation 03 of the ^-distribution as a function of the nominal 
Coulomb-repulsion E c of the fragments at the time of the breakup, i.e. as a function of the 

1 /3 

Coulomb potential of three touching spheres with radii B4 = 1.4 • A { : 



E P . 



E 



Z{ ■ Zj 



1.4 • (A 1/3 + A) 



ll/3-s 



(4) 



This is a generalization of the well known Viola formula |JT|. For events with only two 
heavy fragments the kinetic energy and the Coulomb repulsion are calculated accordingly. 
The experimental results are plotted in figure |9| for the four targets used. Within the 
statistical uncertainties, no target dependence is apparent. In all further plots, mean values 
of the kinetic energy and of the width of the energy distribution for the combined data of 
all four targets will therefore be shown. (E 2 ), (E 3 ), and 02, 03 depend linearly on E c and 
are parameterized in terms of straight line fits (y = m ■ x + h) common to the data of all 
four targets. The slopes and intercepts of these fits are listed in the following table: 





binary 


ternary 


m E 

b E (MeV) 
m,r 

MMeV) 


0.43 ± 0.05 
39.0 ± 4.0 
0.0 ± 0.05 
28.0 ± 3.0 


0.37 ± 0.04 
76.0 ± 5.0 
-0.07 ± 0.01 
44.0 ± 4.0 



The parameters b E and h a describe the mean energies and their variations in the limit 
of E c —0, i.e. without Coulomb interaction, both for binary and ternary events. Under 
the assumption of a purely thermal source with a temperature T and without Coulomb 
interaction the mean values (£'2) and (E 3 ) of the kinetic energy distributions are 2T with a 
width o"2 of a/2T and AT with a width 03 of 2T in case of surface emission and 3/2 • T with 



a width 02 of J 3/2T and 3T with a 03 of y/3T in case of volume emission of the fragments. 
For both breakup scenarios, the temperatures deduced from these relations are within the 
experimental errors identical for binary and ternary events. The assumption of volume 
emission leads to a temperature of 25 MeV whereas the value for surface emission is 20 MeV. 
Results obtained in the reaction Au + Au at 1000 MeV/nucleon where kinetic temperatures 
were extracted from the energy spectra of light charged particles up to 4 He emitted from the 
target spectator [|lt| and temperatures extracted from transverse momentum distributions 
at 600 MeV/u are of similar size ( 15-20 MeV ). 

In line with previous studies |l7"f , the reduced relative velocity is defined as 



Vred(?,3, 



Vrel(hj] 



Z{ + Zj 



(5) 
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where v re i(i, j) is the relative velocity of particles i and j and Zt and Zj are the corresponding 
charges of the fragments. With this definition, the mutual Coulomb repulsion within a 
fragment pair is charge independent. For ternary events, the reduced relative velocity of 
the second and third largest fragment is calculated. Its mean experimental value, averaged 
over all targets, is 0.0206-c ± 0.0005-c. This value will be used later on to adjust the input 
parameters of model calculations. 

The third quantity characterizes the configuration of the three velocity vectors 

«A = ^ , (6) 

where A123 denotes the area of the triangle with its three sides given by the three relative 
velocities v re i(l,2), tv e /(2,3), and v re i(l,3). The normalization A represents the area of an 
equilateral triangle with an circumference of 

u= 1^(1,2)1 + 1^(2,3)1 + 1^(1,3)1 , (7) 

which is the largest area possible for a given circumference. Thus, varies between and 
1, where £7a=0 corresponds to a streched configuration with the three relative velocities 
being collinear and f2A=l corresponding to a situation where the three CM- velocities point 
to the corners of an equilateral triangle. The normalized experimental distributions of the 



reduced area Qa are shown in figure [U] for the four targets: The probability to find an 
equilateral velocity configuration is two orders of magnitude larger than that for a stretched 
one. Within the statistical errors, the distributions are independent of the target, therefore 
the mean value averaged over all four targets was determined to increase the statistics 
especially for small values of In order to address the question of possible correlations 
between the event shape and the charges of the fragments, the average charges (Zi) (z=l,2,3) 
of the three fragments ordered according to their sizes are studied as a function of f^A for 
the combined data of all targets. The results are shown in figure |ll[ Within the statistics, 
the average charges are independent of f^A, indicating that the probability distribution of 
Qa is not a trivial consequence of the charge distribution or the spectator size. 



C. Sensitivity of the three-particle variables 

In order to illustrate the potential sensitivity of the chosen observables, calculations with 
the schematic SOS-code |3f| were performed. This code was especially developped to study 



the influence of two extreme breakup mechanisms on experimentally observable kinetic quan- 
tities, using in both cases a nuclear system of a given size and excitation energy and identical 
multifragment channels. It produces multifragment events with two sets of momentum dis- 
tributions, simulating for each event on the one hand a sequence of binary decays and on the 
other hand a simultaneous breakup using the final partition of the sequential decay chain 
and placing the fragments randomly but without overlap in a sphere. 

For this investigation, masses and excitation energies of the decaying spectator nuclei 
were chosen according to [Q where the authors adjusted the input parameters of a statistical 
fragmentation model (Berlin model) until the relation between (M im f) and Z^^ was well 
reproduced for the system Au + Cu at 600 MeV/nucleon. Since the main motivation 
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of the calculations using the SOS code was to illustrate the potential usefulness of the 
presented observables and not to describe the dynamical aspects of the data, no further 
attempt was done to optimize the input parameters of the code. The standard built-in 
parameters |36[] were used, especially a density for the simultaneous breakup szenario of one 



half of normal nuclear density which is much larger than the values extracted from statistical 
multifragmentation models { p/po = 0.3 in the Copenhagen and the Moscow code and 0.135 
in the MCFRAG code ). 

If the sensitivity of the chosen observables is to be tested, it is - however - important that 
the simulations provide a sample of Monte Carlo data which matches - with respect to the 



fragment composition - the experimental data. This is demonstrated in figure [II], where for 
both breakup scenarios the mean charges (Zi) in ternary events, ordered according to their 
sizes, are compared to the experimental data. The large fluctuations for the simultaneous 
breakup scenario are due ro the fact that only very few events with small values are 
produced ( see next figure ). 

In figure [12|, the probability distribution of the quantity is shown for both breakup 
scenarios and the experimental data. As a reference, the ^-distribution for a thermal sys- 
tem containing three non-interacting fragments is included. For the simultaneous breakup, 
the probability of stretched velocity configurations - i.e. small - is significantly smaller 
than for a purely sequential decay process and for the limit of a thermal system. This dif- 
ference was to be expected, since the repulsive mutual Coulomb interaction shifts initially 
stretched velocity configurations to larger values of The influence of the Coulomb in- 
teraction is especially strong for the relatively small radius used in this simulation, which 
is already an indication that smaller densities will lead to a better description of the ex- 
perimental data. Only due to this repulsion, the velocity configuration is an image of the 
breakup configuration in the coordinate space. Any thermal motion - i.e. any motion which 
is independent of the relative positions of the fragments - reduces this correlation. For re- 
alistic input parameters of the decaying system (see section [TV| ), the correlation coefficient 
r(f2A, ATa) between f^A and the equivalent quantity in the coordinate space Aa 

r(n x ) - {nA • Xa) - {nA) ■ {Xa) (R) 

r( ^ A ' XA) a(n A ) ■ a(X A ) (8) 

has values of approximately 0.1. ( Note that even in the case of T = and three identical 
charges this coefficient does not reach the value 1.0 since the relation between the distance 
of two charged particles and their relative momentum due to the Coulomb-repulsion is 
not linear. ) If - on the other hand - a selfsimilar radial flow dominates the momentum 
distribution, r(f2A, Aa) can reach values around 0.3. 

In figure [13|, the probability distribution of the reduced relative velocity v re d(2, 3) between 
the second and third largest fragments is shown, again both for the data and the SOS 
calculations. The two scenarios predict significantly different relative velocity distributions 
which in both cases differ clearly from the data. In particular, the sequential calculations 
(solid histogram) show a pronounced peak at v rec i(2,3) = 0.012 c. This structure originates 
from the direct splitting of an intermediate state into the observed fragments 2 and 3 at 
a rather late stage of the decay sequence. The absence of this structure in the data may 
therefore signal either a smearing of the relative velocity between the final fragments 2 and 
3 by decays following the splitting into the primordial second and third largest fragments, or 
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proximity effects caused by the presence of other particles, or a different decay mechanism 
which does not produce fragment 2 and 3 via a binary splitting. 

The results presented in figures [L^ and [L3| suggest that the quantities chosen to describe 
the dynamics of the multifragment events are sensitive to important characteristics of the 
decay process. In the following chapter, the experimental results will be compared to cal- 
culations with statistical multifagmentation models and classical three-body calculations in 
order to limit the parameter space of the break-up scenario. 



IV. COMPARISON TO MODEL CALCULATIONS 
A. Statistical multifragmentation models 

Since statistical multifragmentation models have been shown to describe the observables 



in the partition space of the multifragmentation process P,pH^|13l| , it is the obvious next 
step to compare their predictions to the kinetic energy distribution obtained in the present 
experiment. ( It should be emphasized that the description of the partition space comprises 



the cross sections for binary and ternary events as defined in [111 A| . ) Results are shown for 



the Berlin code (MCFRAG) as well as the Copenhagen and the Moscow code. A detailed 
description of the differences between the three models can be found in reference |37| . An 



extensive and detailed investigation of all dynamical observables as defined in sections [III A 
and [llllj| was only performed using the statistical multifragmentation code MCFRAG pfl. 



All three models assume an equilibrated source with a given number of nucleons A at 
a density p with an excitation energy E* per nucleon. This source is non-homogeneous, it 
consists of regions of liquid with normal nuclear density and regions of gas. To compare 
the calculations to the experimentally observed decay of the projectile spectator, the global 
parameters A and E* have to be provided as a function of the impact parameter b. To do 
so, the number of nucleons of the projectile spectator was calculated within a geometrical 
abrasion picture for the collision Au+Au using a radius parameter of 1.3 fm. The excita- 
tion energy for a given spectator size within the three codes was then chosen according to 
|5|-[7|]. For the nuclear density at freeze out the standard values of the models were taken, 
i.e. p/po = 0.3 for the Copenhagen and the Moscow code and 0.135 for the MCFRAG code. 



In figure [14], the size of the projectile spectator and its excitation energy are shown versus 
the impact parameter. It should be noted that for all three models the excitation energy 
necessary to describe the partition space of the multifragmentation is significantly smaller 
than the experimental results obtained for the reaction Au + Au at 600 MeV/nucleon using 
a total energy balance |30|] . The number of events to be produced for each interval in b was 



chosen according to the geometrical cross section for the interval, dP{b) ~ bdb. The impact 
parameter was varied between 0.5 and 12.0 fm in steps of 0.5 fm. For the MCFRAG code, 
the calculation of the observables was done twice: First, the output of the simulations was 
used directly, then random errors on the order of the experimental uncertainties for light 
particles were added to the masses and momenta of the fragments before the same analysis 
was performed. In this way, an upper estimation for the uncertainties produced by the 
experimental resolution was achieved. 

In figure [15] the mean kinetic energies (E2) and (E3) for binary and ternary events ( as 
defined in fll A| ) in the center of mass frame of the two or three particles and the widths of 
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these distributions a% and 03 are plotted versus the nominal Coulomb energy. In the case 
of the MCFRAG code, the results including the experimental resolution are shown, for the 
two other sets of simulations, the uncertainties due to the experimental errors were added 
quadratically to the intrinsic widths of the energy distributions. The mean kinetic energy 
(E3) is reasonably well described by all models, although small differences arise: For the 
whole range of E c , the calculations using the Copenhagen model is steeper than the exper- 
imental distribution, therefore the agreement is, compared to the two other models, worse. 
The overall agreement of data and the three sets of calculations - independent of internal 
details in the theoretical treatment of the fragmentation process - and the simultaneous 
description of (E 2 ) and (E 3 ) by the MCFRAG code are nevertheless a confirmation for the 
expansion of the nuclear matter prior to its decay. The width of the energy distributions, 
on the other hand, is underestimated by almost a factor of two both for events with two 
and three heavy fragments in the exit channel. In spite of deviations between the three sets 
of calculations, the inadequate description of <t 3 is a generic problem of all three statistical 
multifragmentation models. Using the MCFRAG code, it was verified that this underpredic- 
tion of 03 can not be compensated by reasonable fluctuations of the initial excitation energy 
of a given spectator: Combining the events from three sets of calculations with 0.9, 1.0 
and 1.1 times E*(A ) does not change the width of the energy distribution. This variation 
of the excitation energy corresponds within the relevant range of spectator sizes approxi- 
mately to the width of the energy distribution used in || to describe the experimental charge 
distributions. 

In figure [16] the probability distribution for the quantity Qa is plotted both for the data 
and the calculations with the MCFRAG code. The calculated distribution is significantly 
steeper than the experimental one. On the other hand, it is less steep than the result of 
the SOS-calculation for a simultaneous breakup presented in figure |T2|. Since in both cases 
the excitation energy transferred to the spectator matter of a given size is identical and the 
breakup pattern is on average very similar, any differences in the velocity distributions are 
caused by the different radii of the breakup volume. This will result in different contributions 
from the Coulomb interaction and - more important - in different spacial breakup configu- 
rations. On average, an elongated structure will result in a smaller value of than a more 
compact one. If, however, the volume is very small like in the case of the SOS-calculations, 
elongated configurations are less likely. The probability distribution of is therefore ex- 
pected to be steeper than for the more dilute system used for the MCFRAG-calculations. 



B. The influence of angular momentum 

The simulations presented in the previous section showed that the experimental energy 
distributions can not be explained in a purely thermal description of the nuclear matter, 
if the temperature is adjusted to reproduce the charge distributions. It was shown earlier 
that the coupling of random and collective motion increases the fluctuations of the kinetic 
energy |[38| . As an additional degree of freedom angular momentum was therefore taken 
into account. It is well known from the study of fission and compound nuclei at lower 
energies that in heavy ion reactions very large angular momenta can be transferred, causing 
a collective rotation of the excited matter. INC-calculations at 100 and 200 MeV/nucleon 
show that the mean angular momentum per nucleon transferred can be as large as 0.75ft, 
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but even more important than the mean values are the huge angular momentum fluctuations 



which may reach 0.5h per nucleon FWHM [3J|. The influence of angular momentum on the 



decay pattern of nuclear matter within the framework of statistical fragmentation models 
has only barely been studied so far. 

Calculations with the MCFRAG model were done using a version of the code where the 
treatment of angular momentum was implemented in a fully microcanonical way 0. The 
impact parameter was again varied between 3.0 and 12.0 fm in steps of 0.5 fm ( below 5 fm, 
no events with three heavy fragments are produced ) and a total number of 570000 events 
for each set of simulations was produced. In this implementation, the rotational degrees 
of freedom are assumed to be completely thermalized and the contribution of the intrinsic 
rotation of the produced fragments to the total angular momentum is neglected. This is 
supposed to be a good approximation for expanded systems at the time of freeze out, since 
the main part of the angular momentum is contained in the orbital motion of the fragments 
around the common center of mass. 

Calculations were performed for three nuclear densities 0.055p , 0.080p , 0.135p , us- 
ing the relations between impact parameter, system size and excitation energy which were 



already shown in figure [14], and a mean angular momentum (L) of 0.75/iA. The angular 



momentum transfer was distributed according to 

P{L ^ whs - e A^m) ■ (9) 

In reference ||, it was already shown that simulations with this angular momentum distri- 
bution together with a nuclear density of 0.08po describe simultaneously the quantities (E 3 ) 
and CT3. The results, again including the influence of the experimental uncertainties, are 
shown in figure O for the three densities listed above. As expected, the mean kinetic energy 



as well as the width of the energy distribution increases with increasing nuclear density. Due 
to the fact that the mean rotational energy is not very large, the incorporation of angular 



momentum does not change (E 3 ) very much, as a comparison to figure [15] demonstrates, but 
the large variation of angular momenta produces nonthermal fluctuations which increase the 
value of 03 significantly, resulting in a good description of both (E 3 ) and cr 3 for densities 
between 0.055po and 0.080po- At the same time, the quantity Qa is much better described, 
as is shown in figure [18] where the probability distribution of is plotted for the three 
densities. Independent of the nuclear density chosen the probability for the occurence of 
stretched configurations of the three velocity vectors is enhanced. 

To check whether the decay pattern is changed by the angular momentum, observables 
which were used in earlier papers to describe the charge partition space of the reaction 
were investigated: The mean values of the asymmetries 

Z\ — Z2 Z2 — Z 3 

an = „ , „ and a 23 = , (10) 

Zj\ -+- Z/2 ^2 ~T ^3 

the mean number of intermediate mass fragments M im f and the average charge of the 



largest fragment Z max are calculated as a function of Z bound . In figure [19], the results are 
shown for simulations with and without angular momentum together with the experimental 
data. Whereas the mean number of intermediate mass fragments (Mj m /) does not change 
very much under the influence of angular momentum, this is not true for the details of 
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the decay pattern of the spectator: The mean asymmetry (012) between the charges of the 
largest and the second largest fragment decreases dramatically for values of Z bound above 50, 
which means that the two fragments become more comparable in size. As a consequence, 
the mean charge of the largest fragment (Z max ) within an event also decreases. At the same 
time, the mean asymmetry between the charges of the second and the third largest fragment 
(023) increases, which means that in the presence of angular momentum the charge of the 
spectator is more evenly divided between the two largest fragments. The changes in the 
breakup pattern are more pronounced for a small freeze out density. These results are in 
qualitative agreement with the investigations presented by Botvina and Gross ||, where the 
size of the largest fragment and the relative size of the two largest fragments was studied 
under the influence of different amounts of angular momentum. 

From the calculations presented above it is obvious that large angular momenta per 
nucleon destroy the agreement between the results of the statistical multifragmentation 
code and the data as far as the partition pattern of the spectator matter is concerned. This 
is especially true for large values of Z bound , i.e. for peripheral collisions. On the other hand, 
it was shown that the additional degree of freedom increased the fluctuations of the kinetic 
energy by a substantial amount. The question therefore arises whether a better overall 
agreement can be achieved if the transfer of angular momentum per nucleon to the system 
is reduced for large impact parameters. 

If the peripheral reactions are treated in the abrasion-ablation picture applying the for- 
malism described in [fiO|l , values for the angular momentum transfer are obtained which are 
smaller than the value of 0.75/i/nucleon by a factor of 5 to 10. These numbers together with 
a density of 0.135po result in a reasonable description of the partition but the energy fluc- 
tuations are again underestimated. The mean values of the asymmetries ayi and a 2 3 might 
suggest that this can be compensated by an increase of the nuclear density at breakup. 
Unfortunately, this is in contradiction to the description of the quantity (Mj m /). The prob- 
ability to find large values of (M im f) for the Z boun d range between 40 and 70 decreases with 
increasing density. As the mean multiplicity of IMF's is already too small, a further increase 
of the density would make the deviations even worse. 

This leaves no room for a parametrization of angular momentum transfer and density 
which fits both aspects of the experimental data. The charge partition space and the dy- 
namics of multifragmentation events can not be described simultaneously by the statistical 
multifragmentation model even if angular momentum as an additional degree of freedom 
and therefore as a potential source for fluctuations is taken into account. 

The conclusions drawn in this section are valid only for a nuclear system where all 
degrees of freedom are completely thermalized. If this is not true, i.e. if the time scale 
for the equilibration of the rotational degrees of freedom is large compared to that of the 
thermalization of the excitation energy, the process of fragmentation is decoupled from the 
angular momentum transfer. In this case, the amount of angular momentum transferred to 
the spectator does not influence the partition space of the reaction, it only contributes to 
the final momentum distribution of the fragments. Therefore, density and excitation energy 
on the one hand and angular momentum on the other hand can be adjusted independently 
and a reasonable agreement with the experimental data can be achieved. This approach has 
been adopted by the Multics/Miniball group [[41]. It has to be stated, though, that with this 
modification the fragmentation process is not treated in a purely microcanonical picture any 
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longer. 



C. Classical three-body calculations 

A collective radial motion of all constituents of the spectator is another conveivable 
source of fluctuations of the kinetic energy. If the nuclear matter is compressed in the initial 
stage of the reaction, an additional nonequilibrated collective contribution to the motion of 
the nuclear matter will be present |4^j43|| . Even though this effect is expected to be small 
in the peripheral collisions discussed in this paper, values for the radial flow energy up to 
1.5 MeV can not be ruled out ||. First attempts have been made to include collective 
radial flow in statistical models |[44]|, but a consistent implemention is not yet available. 



Therefore, classical three-body calculations were performed to get a quantitative estimate 
for the influence of collective flow. 

The simultaneous emission out of a given volume is modeled in the following way: The 
centers of three non-overlapping fragments with a radius of 1.2- A 1 / 3 are distributed randomly 
within a sphere of radius R. To each fragment, an isotropically distributed initial velocity is 
assigned. Constrained by momentum conservation, these velocities were selected according 
to a probability distribution for the relative kinetic energy 

P(E)~E°.exp(-|) (11) 

with a equals 0.5 or 1.0, corresponding to a volume or a surface emission of the fragments. 
In addition to this random motion, an initial radial flow velocity 



was added to the random velocities of the thermal motion. Here, d± is the position of the 
center of fragment % with respect to the center of mass, e/ is the flow energy per nucleon for 
fragments located at d% = R. The charges and masses of the fragments were obtained by a 
Monte Carlo sampling of the experimental events, thus reducing significantly the uncertain- 
ties associated with the fragment distribution. In order to account for the recoil from light 
particles emitted sequentially from the initial fragments (Z^A^), the measured charges Zi 
and masses Aj were multiplied by a factor (1 — T 2 /(a ■ A)) -1 . For this correction, a level 
density parameter of a=10 MeV was used. The quantity A represents the average energy 
removed by the emission of a nucleon. For simplicity, A = IT + E s + Eb was assumed, where 
E s =8 MeV and -E&=4 MeV are the typical separation energy and barrier height, respec- 
tively. After the interaction of the primordial fragments A^) has ceased, the sequential 
emission of light particles leading to the observed masses and charges [Z^Ai) was assumed 
to take place. For each event, the temperature parameter T was chosen according to the 
experimental value of Zb OU nd from the relation 



T = f T ■ p ■ (79 - Z bound ) , (13) 

where fx is a free parameter. For fj< = 1, the relation describes within the relevant range of 
Zbound reasonably well the temperatures of the initial projectile spectators as predicted by 
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microscopic transport calculations [[16|,[|I|,[| . A value of 0.75 is in agreement with experimen- 
tal results obtained by the He-Li isotope thermometer The paths of the fragments 
were calculated under the influence of their mutual Coulomb field and two-fragment prox- 
imity forces according to Ref. fl47| . Since for the further analysis those trajectories were 



rejected for which the fragments overlapped during the propagation, the influence of the 
proximity force turned out to be rather small. 

In a first step, these schematic trajectory calculations were performed with input pa- 
rameters corresponding on average to those of the statistical model MCFRAG, i.e. a = 0.5, 
fx ~ 0.6 - 0.8, R ~ 7 - 9 fm and €f = 0. The results for (E 3 ) and 03 are comparable 
to those of the statistical model calculations, especially the width 03 is again significantly 
underpredicted in this case. In order to demonstrate this, the schematic calculations for 



fx = 0.7, R = 8 fm and e/ = are included in figure |15|. The agreement of the classical 
calculations and the statistical model calculations for a similar set of external parameters is 
a consistency check and shows in addition that the neglection of the influence of the lighter 
particles produced in the reaction i.e. the restriction of the experimental investigation to 
the two or three heaviest fragments does not change the results significantly. The probabil- 
ity distribution of the quantity is also compared to the results of the statistical model 
calculation ( see figure [TBI ). 



In a next step, the quantities R, ej and fx were varied to fit the experimental data. In 
order to quantify the agreement between the simulations and the experimental observations, 
a reduced x 2 was calculated for each parameter set: 

^ = 1^(5*^0! (14) 

i=l °i 

Here, tu, are the four coefficients characterizing the fits to the three-particle data in figure 
|9] and, in addition, the mean reduced velocity between the two lighter fragments as shown 
in figure [13|. 5i and //j denote the experimental uncertainties of these quantities and the 
corresponding model predictions, respectively. The result is shown in figure ^fj A clear 
minimum of x 2 can be determined for each given flow parameter e/ by varying independently 
the other two model parameters R and fx- The left part of figure ^0] shows in a R - fx 
plane the contour lines with x 2 = 2 for e/ = (R ~ 15 fm), 0.5 (R ~ 22 fm) and 1 MeV 
(R fa 26 fm) and for the two values of the exponent a. The corresponding minima of the 
X 2 -distribution are displayed in the right panel of figure ^ as a function of 6f. Both for 
volume emission and surface emission, values for the flow parameter 6f larger than 1 MeV are 
ruled out whereas the results obtained by values between and 1 MeV show no significant 
difference in Xmin- To demonstrate the quality of the parameter adjustment, the quantities 
(E 3 ), a 3 and v rerf (2, 3) are shown in figure pT| for the parameter set R = 22 fm, ej = 0.5 MeV 
and fx = 1-2. In the lower right part of figure ^l|, f^A - which was not used in the fitting 
procedure - is compared to the experimental values. As expected from the results shown 



in figures 12 and [L6l the probability for the existence of stretched velocity configurations 



increases with increasing radius of the decaying system. The f^A-distribution is nevertheless 



not directly comparable to those shown in figures [12] and [16]: They were achieved assuming a 
fixed breakup density for all decaying systems, whereas the three-body calculations assume 
a fixed breakup volume. 

This set of simulations suggests the disintegration of a highly excited and rather extended 
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nuclear system and very low values of the flow parameter. For Z{, otiri d=55 - the mean value 
for events with three heavy particles in the exit channel of the reaction - the fit values 
correspond to a temperature parameter of approximately 10 MeV and a density below 
0.05po which is much smaller than the values used for the MCFRAG calculations in order 
to reproduce the partition space of the reaction. 

In the framework of these schematic calculations the large freeze out radius is due to the 
balance between Coulomb energy and temperature: If a higher nuclear density is assumed, 
the Coulomb repulsion is much stronger and requires therefore a compensation by a lower 
temperature parameter and a vanishing flow to describe the energy spectra. The fluctuations 
03 of the kinetic energy, on the other hand, reflect in addition to thermal fluctuations 
also fluctuations due to the position sampling within the breakup volume. Thus, lower 
temperatures and especially smaller radii lead to a significant reduction of 03 which can not 
be compensated by the small values of radial flow consistent with the energy spectra. 



V. CONCLUSIONS 

Kinematic correlations between two and three heavy projectile fragments produced in Au 
induced reactions at E/A = 600 MeV have been studied. A comparison of the observables 
to the results of the schematic SOS-model confirms their sensitivity to the disassembly con- 
figuration. Classical trajectory calculations sampling the experimental charge distribution 
limit significantly the possible parameter space of the breakup scenario. Taken at their face 
values these simulations require highly excited and rather extended nuclear systems at the 
time of the breakup. These source parameters differ significantly from breakup parameters 
needed by statistical multifragmentation models in order to describe the observed fragment 
distributions and mean values of the kinetic energy distributions. On the other hand, these 
models are not able to reproduce the fluctuations of the energy distribution. Binary events 
not attributed to binary fission also show fluctuations of the relative kinetic energy which can 
only be described by the same rather high - and probably unrealistic - thermal contribution. 
The introduction of angular momentum into the statistical model improves the description 
of the energy fluctuations, but does not allow anymore to reproduce simultaneously the 
charge partition. 

For any further attempt to reconcile the kinetic observables and the partition pattern of 
the spectator matter two possible approaches seem conceivable: Either the assumption of a 
global equilibrium established prior to the fragmentation process is oversimplified and has 
to be given up, or the statistical models have to be refined. 

The nuclear interaction during the breakup process, for example, is so far ignored, the 
interaction between the fragments is limited to the Coulomb repulsion. ( In the classical 
three-body trajectory calculations present in this work, a nuclear proximity potential is 
included, but its influence is strongly suppressed by the requirement that the fragments do 
not overlap. ) One might speculate that in the case of a stronger overlap of the fragments 
in an earlier stage of the breakup, the nuclear attractive force between the fragments may 
partially compensate the Coulomb repulsion. Thus, smaller radii would not necessarily 
lead to an overestimation of the kinetic energies. First steps to add the nuclear interaction 
between the fragments to statistical decay models in a consistent manner have already been 
undertaken [FISLfSl. A recent publication suggests that the nuclear interaction is indeed 
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relevant for excitation energies up to approximately 10 MeV/u ||50|| . At the same time, 
large fluctuations - similar to dissipative phenomena and shape fluctuations known to be 



important in binary fission [5T| - may arise. 

On the long term, a quantitative understanding of fluctuations and their development 
during the disassembly phase clearly requires dynamical transport models which include a 
realistic treatment of fluctuations on a microscopic level. Significant progress in the devel- 
opment of microscopic transport models has been achieved during the last decade [Q , but 
only recently first microscopic calculations were published which reproduce for the ALADIN 
data both the multiplicity of the fragments and the slopes of their kinetic energy spectra 



53| . In the framework of this model - and in line with previous studies pH- p6|1 - it is found 



that the decaying system is not in thermal equilibrium and that the breakup is dominated 
by dynamical processes. However, the fragment composition agrees with the experimental 
one only for a short time interval after the collision (60 fm/c) and is drastically altered 
during the further time evolution. Thus, a consistent description of the time evolution from 
the first stages of the collision via the formation of primordial excited fragments to their 
eventual deexcitation and formation of individual quantum states within one microscopic 
model is still not available. First attempts to take into account the quantal nature of the 
nuclear system are being pursued |57|-|55] for which the present data may serve as a valuable 
testing ground. 

The authors thank D.H.E. Gross and A.S. Botvina for providing us with their statistical 
multifragmentation codes and for helpful discussions. This work was partly supported by the 
Bundesministerium fur Forschung und Technologie. J. P. and M.B. acknowledge the financial 
support of the Deutsche Forschungsgemeinschaft under the Contract no. Po 256/2-1 and 
Be 1634/1-1. 



21 



REFERENCES 



[1] J. Hubele, P. Kreutz, J.C Adloff, M. Begemann-Blaich, P. Bouissou, G. Imme, I. Iori, 

G. J. Kunde, S. Leray, V. Lindenstruth, Z. Liu, U. Lynen, R.J. Meijer, U. Milkau, 
A. Moroni, W.F.J. Miiller, C. Ngo, C.A. Ogilvie, J. Pochodzalla, G. Raciti, G. Rudolf, 

H. Sann, A. Schuttauf, W. Seidel, L. Stuttge, W. Trautmann, A. Tucholski, Z. Phys. A 
340, 263 (1991). 

[2] P. Kreutz, J.C. Adloff, M. Begemann-Blaich, P. Bouissou, J. Hubele, G. Imme, I. Iori, 
G.J. Kunde, S. Leray, V. Lindenstruth, Z. Liu, U. Lynen, R.J. Meijer, U. Milkau, A. Mo- 
roni, W.F.J. Miiller, C. Ngo, C.A. Ogilvie, J. Pochodzalla, G. Raciti, G. Rudolf, H. Sann, 
A. Schuttauf, W. Seidel, L. Stuttge, W. Trautmann, A. Tucholski, Nucl. Phys. A556, 
672 (1993). 

[3] A. Schuttauf, W.D. Kunze, A. Worner, M. Begemann-Blaich, Th. Blaich, D.R. Bow- 
man, R.J. Charity, A. Cosmo, A. Ferrero, C.K. Gelbke, C. GroB, W.C. Hsi, J. Hubele, 
G. Imme, I. Iori, J. Kempter, P. Kreutz, G.J. Kunde, V. Lindenstruth, M.A. Lisa, 
W.G. Lynch, U. Lynen, M. Mang, T. Mohlenkamp, A. Moroni, W.F.J. Miiller, M. Neu- 
mann, B. Ocker, C.A. Ogilvie, G.F. Peaslee, J. Pochodzalla, G. Raciti, F. Rosenberger, 
Th. Rubehn, H. Sann, C. Schwarz, W. Seidel, V. Serfling, L.G. Sobotka, J. Stroth, 
L. Stuttge, S. Tomasevic, W. Trautmann, A. Trzcinski, M.B. Tsang, A. Tucholski, 
G. Verde, CW. Williams, E. Zude, B. Zwieglinski, Nucl. Phys. A607, 457 (1996). 

[4] C.A. Ogilvie, J.C. Adloff, M. Begemann-Blaich, P. Bouissou, J. Hubele, G. Imme, I. Iori, 
P. Kreutz, G.J. Kunde, S. Leray, V. Lindenstruth, Z. Liu, U. Lynen, R.J. Meijer, 
U. Milkau, W.F.J. Miiller, C. Ngo, J. Pochodzalla, G. Raciti, G. Rudolf, H. Sann, 
A. Schuttauf, W. Seidel, L. Stuttge, W. Trautmann, A. Tucholski, Phys. Rev. Lett. 67, 
1214 (1991). 

[5] Bao-An Li, A.R. DeAngelis, D.H.E. Gross, Phys. Lett. B 303, 225 (1993). 

[6] H.W. Barz, W. Bauer, J. P. Bondorf, A.S. Botvina, R. Donangelo, H. Schulz, K. Sneppen, 

Nucl. Phys. A561, 466 (1993). 
[7] A.S. Botvina and I.N. Mishustin, Phys. Lett. B 294, 23 (1992). 

[8] A.S. Botvina, I.N. Mishustin, M. Begemann-Blaich, J. Hubele, G. Imme, I. Iori, 
P. Kreutz, G.J. Kunde, W.D. Kunze, V. Linddenstruth, U. Lynen, A. Moroni, 
W.F.J. Miiller, C.A. Ogilvie, J. Pochodzalla, G. Raciti, Th. Rubehn, h. Sann, 
A. Schuttauf, W. Seidel, W. Trautmann, A. Worner, Nucl. Phys. A584, 737 (1995). 

[9] A.S. Botvina and D.H.E. Gross, Nucl. Phys. A592, 257 (1995). 
[10] A.S. Botvina, 
A.S. H'inov, I.N. Mishustin, Yad. Fiz. 42, 1127 (1985) (Sov. J. Nucl. Phys. 42, 712 
(1985)). 

[11] J. Bondorf, R. Donangelo, I.N. Mishustin, C.J. Pethick, H. Schulz, K. Sneppen, 
Nucl. Phys. A443, 321 (1985). 

[12] D.H.E. Gross, Zhang Xiao-ze, Xu Shu-yan, Phys. Rev. Lett. 56, 1544 (1986). 

[13] Hongfei Xi, T. Odeh, R. Bassini, M. Begemann-Blaich, A.S. Botvina, S. Fritz, S.J. Gaff, 
C. GroB, G. Imme, I. Iori, U. KleinevoB, G.J. Kunde, W.D. Kunze, U. Lynen, V. Mad- 
dalena, M. Mahi, T. Mohlenkamp, A. Moroni, W.F.J. Miiller, C. Nociforo, B. Ocker, 
F. Petruzzelli, J. Pochodzalla, G. Raciti, G. Riccobene, F.P. Romano, Th. Rubehn, 
A. Saija, M. Schnittker, A. Schuttauf, C. Schwarz, W. Seidel, V. Serfling, C. Sfienti, 



22 



W. Trautmann, A. Trzcinski, G. Verde, A. Worrier, B. Zwieglinski, accepted for publi- 
cation in Z. Phys. A. 

[14] R. Trockel, U. Lynen, J. Pochodzalla, W. Trautmann, N. Brummund, E. Eckert, 
R. Glasow, K.D. Hildenbrand, K.H. Kampert, W.F.J. Muller, D. Pelte, H.J. Rabe, 
H. Sann, R. Santo, H. Stelzer, R. Wada, Phys. Rev. Lett. 59, 2844 (1987). 

[15] D.H.E. Gross, G. Klotz-Engmann, H. Oeschler, Phys. Lett. B 224, 29 (1989). 

[16] Y.D. Kim, R.T. de Souza, D.R. Bowmann, N. Carlin, C.K. Gelbke, W.G. Gong, 
W.G. Lynch, L. Phair, M.B. Tsang, F. Zhu, S. Pratt, Phys. Rev. Lett. 67, 14 (1991). 

[17] Y.D. Kim, R.T. de Souza, D.R. Bowmann, N. Carlin, C.K. Gelbke, W.G. Gong, 
W.G. Lynch, L. Phair, M.B. Tsang, F. Zhu, Phys. Rev. C 45, 338 (1992). 

[18] D.R. Bowman, G.F. Peaslee, N. Carlin, R.T. de Souza, C.K. Gelbke, W.G. Gong, 
Y.D. Kim, M.A. Lisa, W.G. Lynch, L. Phair, M.B. Tsang, C. Williams, N. Colonna, 
K. Hanold, M.A. McMahan, G.J. Wozniak, L.G. Moretto, Phys. Rev. Lett. 70, 3534 
(1993). 

[19] E. Bauge, A. Elmaani, R.A. Lacey, J. Lauret, N.N. Ajitanand, D. Craig, M. Cron- 
qvist, E. Gualtieri, S. Hannuschke, T. Li, B. Llope, T. Reposeur, A. Vander Molen, 
CD. Westfall, J.S. Winfield, J. Yee, S. Yennello, A. Nadasen, R.S. Tickle, E. Norbeck, 
Phys. Rev. Lett. 70, 3705 (1993). 

[20] B. Kampfer, R. Kotte, J. Mosner, W. Neubert, D. Wohlfarth, J.P Alard, Z. Basrak, 
N. Bastid, I.M. Belayev, Th. Blaich, A. Buta, R. Caplar, C. Cerruti, N. Cindro, J.P. Cof- 
fin, P. Dupieux, J. Ero, Z.G. Fan, P. Fintz, Z. Fodor, R. Freifelder, L. Fraysse, S. Frolov, 
A. Gobbi, Y. Grigorian, G. Guillaume, N. Herrmann, K.D. Hildenbrand, S. Holbling, 
A. Houari, S.C. Jeong, M. Jorio, F. Jundt, J. Kecskemeti, P. Koncz, Y. Korchagin, 
M. Kramer, C. Kuhn, I. Legrand, A. Lebedev, C. Maguire, V. Manko, T. Matulewicz, 
G. Mgebrishvili, D. Moisa, G. Montaru, I. Montbel, P. Morel, D. Pelte, M. Petro- 
vici, F. Rami, W. Reisdorf, A. Sadchikov, D. Schiill, Z. Seres, B. Sikora, V. Simion, 
S. Smolyankin, U. Sodan, K. Teh, R. Tezkratt, M. Trzaska, M.A. Vasileiv, P. Wagner, 
J.P. Wessels, T. Wienold, Z. Wilhelmi, A.L. Zhilin, Phys. Rev. C 48, R955 (1993). 

[21] Bao-An Li, D.H.E. Gross, V. Lips, H. Oeschler, Phys. Lett. B 335, 1 (1994). 

[22] O. Schapiro, D.H.E. Gross, Nucl. Phys. A573, 143 (1994). 

[23] P. Glassel, D. v. Harrach, H.J. Specht, L. Grodzins, Z. Phys. A 310, 189 (1983). 

[24] D. Pelte, U. Winkler, M. Biihler, B. Weissmann, A. Gobbi, K.D. Hildenbrand, H. Stelzer, 
R. Novotny, Phys. Rev. C 34, 1673 (1986). 

[25] R. Bougault, J. Colin, F. Delaunay, A. Genoux-Lubain, A. Hajfani, C. Le Brun, 
J.F. Lecolley, M. Louvel J.C. Steckmeyer, Phys. Lett. B 232, 291 (1989). 

[26] G. Bizard, D. Durand, A. Genoux-Lubain, M. Louvel, R. Bougault, R. Brou, H. Doubre, 
Y. El-Masri, H. Fugiwara, K. Hagel, A. Hajfani, F. Hanappe, S. Jeong, CM. Jin, 
S. Kato, J.L. Laville, C. Le Brun, J.F. Lecolley, S. Lee, T. Matsuse, T. Motobayashi, 
J.P. Patry, A. Peghaire, J. Peter, N. Prot, R. Regimbart, F. Saint-Laurent, J.C. Steck- 
meyer, B. Tamain, Phys. Lett. B 276, 413 (1992). 

[27] M. Bruno, M.D'Agostino, M.L. Fiandri, E. Fuschini, L. Manduci, P.F. Mastinu, 
P.M. Milazzo, F.Gramegna, A.M.J. Ferrero, F. Gulminelli, I. Iori, A. Moroni, 
R. Scardaoni, P. Buttazzo, G.V. Margagliotti, G. Vannini, G. Auger, E. Plagnol, 
Nucl. Phys. A576, 138 (1994). 

[28] J. Lauret and R.A. Lacey, Phys. Lett. B 327, 195 (1994). 



23 



[29] G. Bauer, F. Bieser, F.P. Brady, J.C. Chance, W.F. Christie, M. Gilkes, V. Lindenstruth, 
U. Lynen, W.F.J. Miiller, J.L. Romero, H. Sarin, C.E. Tull, P. Warren, Nucl. Instr. and 
Meth. A 386, 249 (1997). 

[30] J. Pochodzalla, T. Mohlenkamp, T. Rubehn, A. Schuttauf, A. Worner, E. Zude, 
M. Begemann-Blaich, Th. Blaich, H. Emling, A. Ferrero, C. Grofi, G. Imme, I. Iori, 
G.J. Kunde, W.D. Kunze, V. Lindenstruth, U. Lynen, A. Moroni, W.F.J. Miiller, 
B. Ocker, G. Raciti, H. Sann, C. Schwarz, W. Seidel, V. Serfling, J. Stroth, W. Traut- 
mann, A. Trzcinski, A. Tucholski, G. Verde, B. Zwieglinski, Phys. Rev. Lett. 75, 1040 
(1995). 

[31] V.E. Viola, T. Sikkeland, Phys. Rev. 130, 2044 (1963). 

[32] R. Brun, F. Bruyant, M. Maire, A.C. McPherson, P. Zanarini, GEANT3 Report, 

CERN/DD/ec/84-1, (1986). 
[33] Th. Stohlker, H. Geissel, H. Folger, C. Kozhuharov, P.H. Mokler, G. Miinzenberg, 

D. Schardt, Th. Schwab, M. Steiner, H. Stelzer, K. Summerer, Nucl. Instr. and Meth. B 
61, 408 (1991). 

[34] T. Rubehn, R. Bassini, M. Begemann-Blaich, Th. Blaich, A. Ferrero, C. GroB, G. Imme, 
I. Iori, G.J. Kunde, W.D. Kunze, V. Lindenstruth, U. Lynen, T. Mohlenkamp, 
L.G. Moretto, W.F.J. Miiller, B. Ocker, J. Pochodzalla, G. Raciti, S. Reito, H. Sann, 
A. Schuttauf, W. Seidel, V. Serfling, W. Trautmann, A. Trzcinski, G. Verde, A. Worner, 

E. Zude, B. Zwieglinski, Phys. Rev. C 53, 3143 (1996). 
[35] G.J. Kunde, PhD thesis (University Frankfurt) 1994. 

[36] J. A. Lopez and J. Randrup, Comp. Phys. Communications 70, 92 (1992). 
[37] D.H.E. Gross and K. Sneppen, Nucl. Phys. A567, 317 (1994). 

[38] U. Milkau, M. Begemann-Blaich, E.-M. Eckert, G. Imme, P. Kreutz, A. Kuhmichel, 
M. Lattuada, U. Lynen, C. Mazur, W.F.J. Miiller, J.B. Natowitz, C. Ngo, J. Pochodza- 
lla, G. Raciti, M. Ribrag, H. Sann, W. Trautmann, R. Trockel, Z. Phys. A 346, 227 
(1993). 

[39] Th. Blaich, M. Begemann-Blaich, M.M. Fowler, J.B. Wilhelmy, H.C. Britt, D.J. Fields, 
L.F. Hansen, M.N. Namboodiri, T.C. Sangster, Z. Fraenkel, Phys. Rev. C 45, 689 
(1992). 

[40] J.-J. Gaimard, K.-H. Schmidt, Nucl. Phys. A531, 709 (1991). 

[41] M. D'Agostino, M. Bruno, N. Colonna, A. Ferrero, M.L. Fiandri, E. Fuschini, 

F. Gramegna, I. Iori, L. Manduci, G.V. Margagliotti, P.F. Mastinu, P.M. Milazzo, 
A. Moroni, F. Petruzelli, R. Rui, G. Vannini, J.D. Dinius, C.K. Gelbke, T. Glasmacher, 
D.O. Handzy, W. Hsi, M. Huang, G.J. Kunde, M.A. Lisa, W.G. Lynch, CP. Montoya, 

G. F. Peaslee, L. Phair, C. Schwarz, M.B. Tsang, C. Williams, A.S. Botvina, P. De- 
sesquelles, I. Mishustin, Proceedings of the XXXV. International Winter Meeting on 
Nuclear Physics, Bormio, ed. I. Iori (Ricerca Scientifica ed Educatione Permanente, 
Milano), 276 (1997). 

[42] J. Hofmann, W. Scheid, W. Greiner, II Nuovo Cimento 33, 343 (1976). 

[43] G. Poggi for the FOPI-Collaboration, Nucl. Phys. A586, 755 (1995). 

[44] Subrata Pal, S.K. Samaddar, J.N. De, Nucl. Phys. A 49, 608 (1996). 

[45] J. Hubele, P. Kreutz, V. Lindenstruth, J.C. Adloff, M. Begemann-Blaich, P. Bouissou, 
G. Imme, I. Iori, G.J. Kunde, S. Leray, Z. Liu, U. Lynen, R.J. Meijer, U. Milkau, 
A. Moroni, W.F.J. Miiller, C. Ngo, C.A. Ogilvie, J. Pochodzalla, G. Raciti, G. Rudolf, 



24 



H. Sann, A. Schuttauf, W. Seidel, L. Stuttge, W. Trautmann, A. Tucholski, R. Heck, 

A.R. DeAngelis, D.H.E. Gross, H.R. Jaqaman, H.W. Barz, H. Schulz, W.A. Friedman, 

R.J. Charity, Phys. Rev. C 46, R1577 (1992). 
[46] W. Bauer, Phys. Rev. C 38, 1297 (1988). 
[47] J. A. Lopez and J. Randrup, Nucl. Phys. A503, 183 (1989). 
[48] L. Satpathy, M. Mishra, A. Das, M. Satpathy, Phys. Lett. B 237, 181 (1990). 
[49] A. Das, M. Mishra, M. Satpathy, L. Satpathy, J. Phys. G 19, 319 (1993). 
[50] C.B. Das, A. Das, M. Satpathy, L. Satpathy, Phys. Rev. C 56, 1444 (1997). 
[51] For a recent review see F. Gonnenwein, The Nuclear Fission Process, Edt. C. Wagemans, 

CRC Press, Boca Raton, 287 (1991). 
[52] See e.g. The Nuclear Equation of State, Part A and B, Edt. W. Greiner and H. Stocker, 

Plenum Press, New York (1989). 
[53] PB. Gossiaux, R. Puri, Ch. Hartnack, J. Aichelin, Nucl. Phys. A619, 379 (1997). 
[54] D.H. Boal, J.N. Gosli, C. Wicentowich, Phys. Rev. Lett. 62, 737 (1989). 
[55] D.H. Boal, J.N. Gosli, C. Wicentowich, Phys. Rev. C 40, 601 (1989). 
[56] G.J. Kunde, J. Pochodzalla, J. Aichelin, E. Berdermann, B. Bethier, C. Cerruti, 

C.K. Gelbke, J. Hubele, P. Kreutz, S. Leray, R. Lucas, U. Lynen, U. Milkau, 

W.F.J. Muller, C. Ngo, CH. Pinkenburg, G. Raciti, H. Sann, W. Trautmann, 

Phys. Lett. B 272, 202 (1991). 
[57] A. Ohnishi and J. Randrup, Phys. Lett. B 394, 260 (1997). 
[58] J. Schnack and H. Feldmeier, Prog. Part. Nucl. Phys. 39, 393 (1997). 
[59] J. Schnack and H. Feldmeier, Phys. Lett. B 409, 6 (1997). 



25 



FIGURES 




FIG. 1. Schematic view of the experimental setup in the bending plane of the magnet. The 
beam enters from the left and hits the beam counters before reaching the target. Mid-rapidity 
particles are detected in the Si-Csl array. Projectile fragments are tracked and identified in the 
TP-MUSIC II detector and in the time-of-flight wall. 
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FIG. 2. Schematic illustration of the design and the operation of the MUSIC II detector. In 
the first and in the last cage, the electric field is orientated horizontally, in the second field cage 
vertically, providing the horizontal, the vertical and again the horizontal track information. After 
amplification and pulse shaping the signals are recorded by a sampling ADC. The insert shows a 
raw time spectrum as delivered by the 16 anodes of each field cage for a track with a finite angle 
with respect to the anode plane. 
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FIG. 3. Charge spectrum for particles detected by the MUSIC detector in the reaction Au+Cu. 
Single charge resolution is obtained for the whole range of fragments from the lower detection 
threshold Z=8 up to beam particles. Fragments with charge 70 and above are suppressed by the 
trigger conditions. 
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FIG. 4. Rigidity R of gold projectiles at 600 MeV/nucleon versus the horizontal position in the 
reference plane of the MUSIC detector. Only events were selected where the gold nuclei passed 
through the carbon target without any nuclear interaction. The three charge states 79 + , 78 + and 
77 + correspond to the rigidities 2.990, 3.028 and 3.068 GeV/c, respectively, which are different by 
1.3% (dashed lines). 
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FIG. 5. Mass spectrum of particles detected in the MUSIC detector for the systems Au+Al 
and Au+Cu. Single mass resolution is obtained for charges up to 12. 
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FIG. 6. (top): Differential cross sections da / dZb OU nd both for the inclusive data (histograms) 
and for events with two and three fragments of charge >8 in the exit channel. For the system 
Au+C the cross section for binary fission is also shown (solid crosses). 

(bottom): Fraction of Zt, oun d contained in the sum of the charges of the two or three heavy fragments 
versus Z^ oun d. The symbols represent a cross section weighted mean value for all four targets. 
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FIG. 7. Mean multiplicity (M/ p ) of light particles from the interaction zone of the reaction 
detected in the hodoscope versus Zb OU nd- The histograms show the inclusive distributions for the 
four targets. 

(top): Binary events where binary fission - as defined in the text - was excluded, 
(bottom): Ternary events 

(Mi p ) decreases with increasing size of the projectile spectator. Within the experimental errors, the 
multiplicity for a given value of Z^ound is independent of the number of the projectile fragments. 
This holds for the whole range of Z^ound with the exception of very peripheral reactions with 
Zbound <65 where the inclusive distributions are dominated by spallation reactions. 
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FIG. 8. Mean transversal velocity fitrans = \J P% + Py of the center of mass relative to the 
velocity of the beam for events with two (upper panel) and three (lower panel) heavy particles 
with Zi >8 in the exit channel. The bounce decreases with decreasing centrality of the reaction. 
The non-zero value for the largest Zf, oun d bins is due to the trigger condition which requests the 
detection of at least one light particle in the hodosope. The dashed lines in both panels show the 
same linear fit to the binary data, to demonstrate that for a given value of Z^und the distributions 
do not depend on the multiplicity of heavy particles. 
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FIG. 9. (left panels): Mean value (£3) (top) and the standard deviation 03 (bottom) of the 
total kinetic energy as a function of the nominal Coulomb energy E c for events with three large 
fragments in the exit channel of the reaction. 

(right panels): (£2) (top) and the standard deviation er 2 (bottom), the equivalent quantities for 
events with two large fragments where binary fission is excluded. 

The results are shown for the systems Au+C, Al, Cu and Pb at E/A = 600 MeV. Both the 
mean kinetic energy and the width of the energy distribution are within the experimental errors 
independent of the target. The straight lines are least square fits to the combined data of all 
targets. 




FIG. 10. Experimental probability distributions for The symbols for the combined data of 
all targets are shifted by 0.02 in x-direction. The event shape in velocity space for large and small 
values of are indicated schematically in the upper left and right corners. 
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FIG. 11. Mean values of Z{ (i=l,2,3) as a function of SIa for the combined data of all targets 
(symbols). Also included are the results of calculations with the SOS-model assuming a sequential 
(solid histograms) and a simultaneous (dashed histograms) decay of the projectile spectator ( see 
section [III C ). For all three fragments, the mean charges are independent of f2 a , showing that the 
probability distribution for Qa is not a trivial consequence of the charge distributions. 
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FIG. 12. Probability distribution of the observable ^a- The stars refer to the combined data 
of all targets. Predictions of the SOS-code assuming a sequential or prompt break-up are shown by 
the solid and dashed histograms, respectively. Also included are the values for a thermal system 
without Coulomb interaction, (filled circles). 
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FIG. 13. Probability distribution of the reduced relative velocity f rec z(2, 3) between the second 
and third largest fragment. The symbols refer to the combined data of all targets. Predictions 
of the SOS-code assuming a sequential or prompt break-up are shown by the solid and dashed 
histograms, respectively. 
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FIG. 14. Input parameters for the simulations with statistical multifragmentation codes, 
(top): Size of the decaying spectator Aq versus the impact parameter b. 

(bottom): Excitation energy E* per nucleon versus the impact parameter. The short dashed, long 
dashed and solid lines show the values used for the three codes COPENHAGEN, MOSCOW and 
BERLIN, respectively. 
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FIG. 15. Mean kinetic energies (E3) and (E2) (top) and the standard deviations u 3 and 02 
(bottom) as a function of the nominal Coulomb-energy E c . The symbols denote the experimental 
data averaged over the reactions Au+C, Al, Cu and Pb. The short dashed, long dashed and solid 
histograms present predictions of the COPENHAGEN, MOSCOW and MCFRAG fragmentation 
models, respectively. The dashed-dotted line shows the result of a classical trajectory calculation 



as described in section IV C . All calculations underpredict significantly the width of the energy 
distribution while the mean value of the kinetic energy is well described. 




FIG. 16. Probability distribution for fi^- The symbols denote the experimental data, the 
histograms model calculations: The solid line represents the calculation with the MCFRAG code 
using the standard input parameters of the model. The dashed-dotted line shows the result of a 
classical trajectory calculation as described in section fVC, where the same input parameters were 
used. 
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FIG. 17. Mean kinetic energy (E3) and the standard deviation 0-3 as a function of the nominal 
Coulomb energy E c under the assumption of a mean angular momentum transfer to the spectator of 
0.75?i/nucleon. The dotted, dashed and solid histograms present MCFRAG-calculations with freeze 
out densities of 0.055po> 0.080/9o and 0.135/9o, respectively. The symbols represent the experimental 
data. 




FIG. 18. Probability distribution for Qa under the assumption of a mean angular momentum 
transfer to the spectator of 0.75/i/nucleon. The dotted, dashed and solid histograms present 
MCFRAG-calculations with freeze out densities of 0.055po 5 0.080/9o and 0.135/?o, respectively. The 
symbols denote the experimental data. 
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FIG. 19. Mean charge asymmetry between the two largest fragments (012) (top left) and be- 
tween the second and third largest fragment (023) (top right), mean number of intermediate mass 
fragments (Mj m /) (bottom left) and mean charge of the largest fragment (Z max ) (bottom right) 
versus Zf, oun( i under the assumption of a mean angular momentum transfer to the spectator of 
0.75?i/nucleon. The stars denote the experimental data. The dotted, dashed and solid histograms 
present calculations with freeze out densities of 0.055po> 0.080/3o and 0.135po> respectively. For com- 
parison, the results of calculations without angluar momentum transfer and a freeze out density of 
0.135/3Q are included as open circles. 
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FIG. 20. Parameter adjustment for the classical three-body calculations. The solid lines corre- 
spond to surface emission (a = 1.0), whereas the dashed lines show the results for volume emission 
(q = 0.5). 

(left): Contour lines for x 2 =2 in a plane defined by the volume radius R and the scaling factor fx 
of the temperature for flow parameters tf = 0.0, 0.5 and 1.0 MeV. 

(right): Minima of the x 2 -distribution as a function of e/. Values of e/ larger than 1 MeV are ruled 
out whereas values between and 1 MeV show no significant differences in 
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FIG. 21. Mean kinetic energy (£3) and its standard deviation 0-3, v re d(2, 3) and for classical 
three-body calculations (histograms). The parameter set R=22 fm, e/= 0.5 and /t=1.2, corre- 
sponding to the minimum in x 2 in the case of surface emission of the three fragments ( a = 1.0 ) 
was chosen. The symbols represent the experimental data. 
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